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Abstract — A clear understanding the behavior of the error 
probability (EP) as a function of signal-to-noise ratio (SNR) 
and other system parameters is fundamental for assessing the 
design of digital wireless communication systems. We propose 
an analytical framework based on the log-concavity property of 
the EP which we prove for a wide family of multidimensional 
modulation formats in the presence of Gaussian disturbances 
and fading. Based on this property, we construct a class of local 
bounds for the EP that improve known generic bounds in a given 
region of the SNR and are invertible, as well as easily tractable 
for further analysis. This concept is motivated by the fact that 
communication systems often operate with performance in a 
certain region of interest (ROI) and, thus, it may be advantageous 
to have tighter bounds within this region instead of generic 
bounds valid for all SNRs. We present a possible application 
of these local bounds, but their relevance is beyond the example 
made in this paper. 

Index Terms — Error statistics, fading channels, local bounds, 
log-concavity, performance evaluation, probability. 



I. Introduction 

THE performance evaluation for digital wireless commu- 
nication systems in terms of bit error probability (BEP) 
and symbol error probability (SEP) requires a careful charac- 
terization of disturbances, such as noise and interference, as 
well as of the wireless channel impairments due to small-scale 
and large-scale fading (see, e.g., [l]-[3]). This can result in 
cumbersome expressions for the error probability (EP) which 
require numerical evaluation^ 

At a first thought, this fact does not appear a relevant 
issue from the performance study point of view due to the 
increasing trend of computational power of computers. On 
the other hand, these cumbersome solutions do not provide 
a clear understanding of the performance sensitivity to system 
parameters, which is of great importance for system design, as 
well as they are often too complicated for further evaluation 
or inversion, which is as example needed in order to obtain 
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' Hereafter, when EP is indicated without specification of BEP and SEP it 
means that the concept is valid for both BEP and SEP. 



thresholds in adaptive communication systems (see, e.g., [4]- 
[6]). Moreover, it has to be emphasized that simple parametric 
approximations and bounds on the performance at lower lay- 
ers, such as physical layer, can avoid long bit-level simulations 
in upper level protocols network simulators, provided that they 
are able to capture the main aspects affecting the performance 
at lower levels. 

Mainly, but not only for these reasons, the derivation of 
approximations and bounds on the exact EP is still of interest 
in the communication theory community. An example is given 
by Af-ary quadrature amplitude modulation (AI-QAM), that 
is adopted in several standards for wireless communication 
systems, due to its bandwidth efficiency, and is largely stud- 
ied in conjunction with adaptive techniques which change 
modulation parameters to maximize transmission rate for a 
given target BEP in wireless channels. In fact, although early 
work on 7\f-QAM dates back to the early sixties [7]-[10], the 
evaluation of BEP for arbitrary M is still of current interest0 
To briefly summarize some relevant results for additive white 
Gaussian noise (AWGN) channel we recall that: parameterized 
exponential approximations fitting simulative BEP are adopted 
in [12]-[15]; approximations based on signal-space concepts 
were given in [16]; an exact method to derive the SEP 
was proposed in [17]; a recursive algorithm exploiting the 
relationship among different constellation sizes was developed 
in [18]; exact expression of the BEP for general M was 
derived in [19]. Comparisons among approximations, bounds 
and the exact solution in fading channels (with small-scale 
fading and large-scale fading, i.e., shadowing) are given in 
[5], where it is shown that, for low and medium values of the 
signal-to-noise-ratio (SNR), approximations depart from exact 
solutions as the constellation size, M, increases. Moreover, 
small differences between exact solution and approximation in 
AWGN channel can become relevant when the instantaneous 
BEP is averaged over small-scale fading. Similarly, in systems 
with multichannel reception, known approximations depart 
from the exact EP as the diversity order increases [20], [21]. 

It is well known that bounds carry more information than 
approximations and also enable system design based on the 
worst or best case. Quite often bounds are tight to the exact 
EP only for high SNRs (namely asymptotic bounds). Here 
we are interested in deriving simple invertible bounds tight 
in a given region of interest (ROI) for the EP (e.g., for the 
BEP of uncoded systems typical ROI's are [10"'^, 10~^] or 

^For a brief history of J\/-QAM, see [11]. 
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[10-4,10-2]). 

In this paper, we define the concept of locally-valid bounds 
(in the following called local bounds), that are tight upper and 
lower bounds on the EP valid only within a given region of 
the EP and not for all SNRs. This concept is motivated by 
the fact that there is often a ROI for the performance of the 
system under consideration and it is preferable to have tight 
bounds in this region instead of bounds valid for all SNRs 
which are far from the exact solution within the ROI. 

The behavior of the EP is important for the definition 
of local bounds. In fact, the proposed framework is based 
on its log-concavity property. We recall that a function z( ) 
is log-concave if logz(-) is concaved In most cases the 
EP is reported on log-scale and investigated as a function 
of the signal-to-total disturbance ratio in deciBel (dB). It 
is commonly recognized that on this scale the function is 
concave in several cases of interest. Even though it is generally 
assumed, as the authors often acknowledge, there is not known 
formal proof of the long-concavity of the EP (examples of 
related issues are: convexity properties in binary detection 
problems which were analyzed in [22], and some results for 
the asymptotic behavior of bounds that were investigated in 
[23]). 

In this paper we introduce the problem of log-concavity 
for general multidimensional decision regions and we prove 
this property for a class of signals with constellation on a 
multidimensional regular grid in the presence of Gaussian dis- 
tributed disturbances, such as thermal noise and interference. 
In fact, there are several wireless systems and situations in 
which the interference can be modeled as Gaussian distributed 
(see, e.g., [24]-[29]). After having proved the log-concavity 
in both AWGN and fading plus AWGN channels for single 
and multiple channels reception schemes, as examples of 
application, it will be shown how to take advantage of this 
property in order to simplify the derivation of bounds valid 
for all SNRs, as well as to define local upper bound (LUB) 
and local lower bound (LLB) valid in a given ROI. Moreover, 
the form of the local bounds and the fact that they are easily 
invertible also enables the derivation of local bounds for others 
relevant performance figures such as the EP-based outage 
probability, which is the probability that the EP averaged 
over small-scale fading exceed a given tolerable target value 
[30], also exploited for the evaluation of the mean spectral 
efficiency for adaptive modulation techniques [6]. It is finally 
emphasized that the log-concavity property for the EP can 
have many others applications, thus its relevance is beyond 
the applications illustrated in the paper. 

The rest of the paper is organized as follows: in Sec. Ullthe 
log-concavity property of the EP is proved in AWGN and in 
AWGN plus fading for systems employing single and multiple 
channels reception, and in Sec.|III]it is applied to define a new 
class of bounds and local bounds, with a discussion on possible 
applications. Finally, our conclusions are reported in Sec. HV] 



^In the paper notation log stands for natural logaiithm. 



II. Log-Concavity Property for the Error 
Probability 

In this section the log-concavity property for the EP is 
discussed first for transmission in AWGN channel by high- 
lighting the mathematical structure of the problem in the 
different applications of digital communications. Since we are 
interested in obtaining general results we base our frame- 
work on the origin of detection errors in the presence of 
Gaussian disturbances. Within this framework we will then 
prove the log-concavity property for the class of signals with 
constellation on a multidimensional regular grid (e.g., in the 
two-dimensional case this class includes the well known M- 
QAM constellation). Finally, we will address the log-concavity 
property in systems with AWGN plus fading channels. 

Consider a set X of M constellation points on a d- 
dimensional signal space, i.e. X = € : i — 

1 , . . . , Af I . Let us consider an arbitrary probability distribu- 
tion on the set X and let pi denote the probability of a point 
Xi E X for 1 < i < M (we arbitrarily order points in X 
without loss in generality). For each Xi G X let us define its 
neighborhood 

Ri = ^x eR"^ : \x - Xi\ = iniii\x ~ Xj\^ (1) 

as the set of points closest to Xi in X. Suppose that we transmit 
a point Xi with probability pi in an AWGN channel, hence 
we receive xi + ag where cr > 0, and g has a standard 
Gaussian distribution 1) on M'* with mean zero and 

identity covariance matrixO We classify each point according 
to a region Ri that it belongs to, which means that we make 
an error if Xi + ag ^ Ri or 

g ^ a-\R, - X,) = {a-\y ~ x,) : y e R,} . (2) 

Through the change of variabl^ a^^ = e*, the total probabil- 
ity of making an error results in @ 

p{t,X) = J2Pi^{9<^e\R,-x,)} 

i 

= l-^KPlgGe^li?, -a;,)}- (3) 

i 

If we denote by /i a standard Gaussian measure on M then 
the distribution of vector g is a product measure /i®'^ and, 
therefore, 

p {t, X) = Y^p, f^'^" (e*i?[) = 1 ~ (e*i?0 , (4) 

i i 

where we denoted by i?'; = Ri — Xi the region Ri translated 
by Xi, and by R^ the region M'' — i?^. 

The function p {t, X) is the error probability in the detection 
of digital signals, either coded or uncoded, as a function of 
signal-to-noise ratio t (in logarithmic scale). Proving the log- 
concavity property of this function is a challenging task. In 

'^Note that cr^ is inversely proportional to the SNR. 

^Which is strictly related to the transformation of the SNR from the linear 
to the dB scale. 

^Notation P {A} stands for probability of event A. 
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k=0 
d 

= l-Y.Pk M([-(«/2)e*,oo))'~V(h(a/2)e*,(a/2)e*])', 



fc=0 



(9) 



fact, the log-concavity of single terms in the linear combina- 
tion of eq. (|4| depends on the specific structure of regions 
R[ and R[ and in any case a possible linear combination of 
log-concave functions is not necessarily log-concave. 

Only in few special cases, as example when all the regions 
i?J have the same measure and a special symmetry around 
the axis intersecting points and ~Xi, these issues may be 
overcome with the help of the Prekopa-Leindler theorems [31], 
[32] which states that the function F{x) = J^f{x.,y)dy, 
where x e M",y e K™, is log-concave in R" ;/ /(x, y) 
is log-concave in and A is a convex subset o/M™. 

Two examples, one for uncoded system and the other for 
coded system, will illustrate these simple cases below. In all 
the other cases an inspection of log-concavity property should 
be based on the specific properties of the signal set X. In 
the next Sec. III-AI we will provide the proof of log-concavity 
property for the specific case of signals with constellation on 
a multidimensional regular grid, which covers all the relevant 
applications based on A/-QAM signaling. 

Example 1 ( M-PSK): Let us consider a 2-dimensional signal 
set X where the M points are regularly placed on a circle. The 
angular separation between closest points is 2tt/M (see Fig. 
[TJleft). This is the signal set used by M-PSK signaling. The 
regions R[ are circular sectors, have the same form and the 
same measure, and are convex. The same holds for regions 
i?';, which are concave instead. If we split each region R'- in 
two parts, and p^^\ using the line connecting the points 
(0,0) and —Xi, all these sub-regions have the same Gaussian 
measure and are convex. Since p(t,X) = 2/1*^^ ^e^r^V 
the log-concavity of a single term has to be checked. The 
2-dimensional Gaussian measure can be evaluated by using 



polar coordinate^ in 



as 



(1) 



1 

27r 



it/2 



(5) 



,(1) 



Since the function e ^^'^ is log-concave 
X [— 7r/2 + tt/A/, 7r/2], the Prekopa-Leindler 



'TT/2+Tr/M 

where S{d) — sui{-k / AI) / cos{6) describes the boundary 
of region p, 
for (t, e) e ' 

theorerr|j assures that p {t, X) is log-concave with respect to 
t. 

Example 2 (Parity check linear block codes and BPSK): Let 
us consider the k + 1-dimensional signal set X representing 
signals obtained by combining a simple parity-check {k + l,k) 
binary block code and binary antipodal modulation. All the 
M = 2*^ points are placed on (half of) the vertexes of a fc + 1- 
dimensional cube and are equidistant from the origin. Each 

'The origin is tlie point Xi and 9 is the angle with respect to the line 
orthogonal to p'^' boundary. 

^Here, the domain is restricted to R X [— 7r/2, n/2]. 





Fig. 1. Decision regions Ri for the examples 1 (left) and 2 (right) of Sec. 
im respectively. 



point has ric = fc(fc+l)/2 closest points or neighbors and each 
region Ri is bounded by faces in the k + 1-dimensional 
space (see Fig. [TJright). All regions R[ have the same form 
and the same measure, and are convex. The same holds for 
regions R'^, which are concave instead. Let us now simplify 
the example to fc = 2 for better understanding. We have 4 
equidistant points placed on 4 vertices of a cube. Regions 
R[ and R'^ are bounded by rtc = 3 planes intersecting in 
—Xi. If we split each region R'- into three parts, p'^^\ pf'^ and 

(3) 

Pi , using three half-planes generated by the line connecting 
the points (0,0,0) and —Xi, all these sub-regions have the 
same Gaussian measure and are convex. Since p {t, X) = 
^eVl^'^ J the log-concavity of a single term has to 

be checked. By using cylindrical coordinate^ in R'^, the 3- 
dimensional Gaussian measure can be evaluated as 

^03 / t (i)\ ^ i Q U(l + SHe))eH/3) de (6) 

where zS{Q) = z j {^/Zcos(e)) with z > ^2/3 de- 
scribes the boundary of region p\^\ (3 = arcsin(-\/2/3). 



and Q{.) is the Gaussian Q-function. Since the function 

Q (v'(l + S'2(6'))e*4/3) is log-concavf^] for {1,9) e R x 



[— 7r/2 + arcsin(^2/3), 7r/2], the Prekopa-Leindler theorem 
assures that p (t, X) is log-concave with respect to t. 

A. Log-Concavity property: signals with constellation on a 
multidimensional grid 

Given a > 0, consider a set X of M — n'^ points on R'' 



X = 



|(fcia, . . .,kda) : h 



= 1, 



for ain < d 



} (7) 



'Here, the origin is the point Xi, z is the coordinate along the line 
orthogonal to boundary, is the angle on the plane orthogonal to z- 
axis. 

'"Note that Q{y/x) is log- concave, whereas (1 -f S^{e))e^ IS convex. 
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PMi) 



'M log2 M 



log2 \/M (1-2"'=)v^-1 

E E (- 

k=l 1=0 



2fc-i _ 



i ■ 2 



erfc 



(2z + l). 



' 37 
2(M- 1) 



(21) 



that form a regular finite grid on with each coordinate 
taking n possible values a, 2a, . . . , na^^ Since A" is a regular 
grid, all sets R[ take a particularly simple form, namely, each 
such set is equal to one of the d + 1 sets given by 

[-a/2,ooY-^ x[-a/2,a/2f for fc = 0, . . . , (8) 

up to a permutation of coordinates. The product measure /.t®'' 
is invariant under permutation of coordinates, thus we can 
identify each set R[ with one of the sets in dHJ. If Pk is the 
sum of probabiUties pi of points contained in the regions of 
type fc, then we obtain (|9]l. 

We note that making a change of variables t — > 
{t — log(a/2)) it suffices to consider the case of a/2 — 1. 
Let us now define 



g := g(e*)^M([-e*,oo)) , 
h:= h{e')^^,{[~e\e']) . 

This leads to the following representation 



p{t,X) = 1 



d 

k=0 



(10a) 
(10b) 



(11) 



Since g = (h + l)/2 we obtain 

p{t,X) = l-H{h{e^)) ^p{n,d,t) , (12) 

where 



k=0 



1 



d-k d 

h'' =Y,PkHk{h). (13) 

fe=0 



Remark. Notice that all derivatives of H with respect to h are 
nonnegative and H{1) = 1. 

We now prove the log-concavity of p {n, d, t) with respect 
to t starting with 2 Lemmas. The main Theorem with proof 
will follow. 

Lemma 1: For any h E [0, 1] and with H{h) given by ( fT3] l 
the following inequality holds 

(1 - h) H' {hf-{l - H {h)) {H' {h) - (1 - /i) H" (h)) > . 

(14) 

Lemma 2: If d > 1, for any h G [0,1] and with Hk{h) 
given by (fTJI i the following inequality holds 

{l-h)H',{h)H'^{h) 

- [1 - Hk{h)] [H'^{h) ^{\-h) H';^{h)] > . (15) 
We will prove this Lemma in the Appendix and now show 
how Lemma [T] implies the main Theorem [T] 

Theorem 1: For any n, d > 1 the function t ^ p{n, d, t) is 
log-concave. 

"without loss of generality X can be translated. 



Proof: (ofTheoremU} Let G (e*) = H {h (e*)). A simple 
calculation gives 



^log(l-G(e*))=-e* 



(16) 



(1 - G (e*)) G' (e*) + e*G' (e*)' + e* (1 ~ G (e*)) G" (e*) 
(1-G(e*)f 
The right hand side is negative if an only if 

(1 - G (c)) G' (c) + cG' {cf + c (1 - G (c)) G" (c) > , 



where c = e* > 0. Since 

G' (c) = i?' (/i) /i' (c) , 

G" (c) if" (h) h' {cf - H' (h) ■ ch' (c) , 

and by definition of h{c) giving 

^2 



(17) 



h" (c) = -c h' (c) , 
we can rewrite ( flTl i as 

(1 - i/) (i/'/i' (1 - c2) + c (/i')' H"^ + c (/i')' (^')' > , 
or, equivalently. 



h' > (c^ - 1) 



\2 ■ 



(18) 



H" {1-H) + {H' 

It is immediate to see that if c < 1 then the inequality holds 
For c > 1 the proof follows from 

[l-H) H' 



>l-h> 



\2 ■ 



- 1 H" {1-H) + {H') 

Here, the left hand side inequality can be derived from 

c 



- 1 



-er— > 



e ^ dx , 



(19) 



(20) 



which is verified for c > 1, and the right hand side 
inequality follows from Lemma 1. □ 

Remark. The instantaneous BEP expression for coherent 
single reception M-QAM systems with arbitrary M as a 
function of the instantaneous symbol SNR 7 is given by 
(I2TI 1 where [xj denotes the integer part of x [19]. One 
might consider to try to prove log-concavity directly using 
this explicit expression. However, this seems to be a difficult 
task since sum of log-concave functions is not log-concave in 

Both sides of )20t tend to as c — > 00, therefore, it is enough and 
simple to show that for any c > 1 



dc V — 1 J dc 



1: 



e 2 dx 
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general and the BEP is a linear combination of positive and 
negative terms containing the complementary error functiorQ, 
erfc(-), making the analysis of (I2TI 1 not at all straightforward. 

B. Log-Concavity property: systems with AWGN plus Fading 

In the above proof of log-concavity of the function t 
p{t, X) the size of the grid a > was fixed. When we transmit 
a symbol related to the constellation point Xi with fixed a in 
AWGN plus fading channel, the receiver observes Fxi + ag 
where F is a random variable (RV) representing the channel 
gain due to fading. This is equivalent to the observation of 
Xi+ag when the constellation has a random scaling parameter 
a with the same statistics of the fading gain Thus, we 
show now that when the constellation X is scaled with a real 
random parameter a, the average of p{t, aX) over a is still 
log-concave in t. 

Let us now make the dependence of p on a explicit, through 
the change of variable b = 2 log a, thus a/ a = e*+''/^. We 
obtain the function p{t + &/2, X), which is log-concave as a 
function of both variables (t, 6) if p{t, X) is log-concave with 
respect to t. Hence, what follows is valid for all signal sets 
X with log-concave instantaneous EP function of t. To obtain 
the EP averaged over fading we have to evaluate the expected 
value of p{t + 6/2, X) with respect to the RV h. 

Theorem 2: If 6 has a log-concave PDF then the average 
of p over 6, that is X) = Ef,{p(i + 6/2, X)}, is also log- 
concave. 

Proof: {of Theorem |2]i Suppose that 6 has a distribution 
with log-concave density, that is density of the form e^^^'') 
for some convex function ^(6). Then'^l 

p{t)^ j p{t + b/1)e-^^''^db 

is the average of p over 6. Since p{t + 6/2)e~^(^) is log- 
concave in both variables (t, 6), Prekopa-Leindler inequality 
[31], [32] impHes that p{t) is log-concave. □ 
Theorem |2] shows that if the distribution of 6 has log- 
concave density then the average over a is log-concave. This 
apply to several cases of interest (e.g., single and multiple 
channel reception in Rayleigh, Nakagami-m, and log-normal 
fading) thus leading to log-concave average EP. This can be 
verified by considering that if the PDF fa^{£,) of is given 
then the PDF of 6 = log results 



(22) 



For Nakagami-m fading (having m > 1/2) the PDF of is 
given b\l'^ 



r(m) 



(23) 



'^The complementary error function is in well known relationship with he 
Gaussian Q-function, i.e., Q{x) = (1/2) erfc |^2:/\/2j- 

'''E.g., this represents the case of flat fading channel and coherent reception. 
'^We omit here the dependence on set X. 

'*lt is well known that Rayleigh fading is included in Nakagami-m when 
m = 1. 



from which by (l22T i we obtain 

h{z) 



T{m) 



(24) 



that is log-concave in z since m{z — e^) is concave and 
m™/r(TO) > 0. For log-normal fading (i.e., in dB is a 
zero-mean Gaussian RV with variance cr^g) the PDF of is 
given by {ly = 10/ log 10) 



from which by 



2vra-dBf 
we obtain 



(25) 



(26) 



'2'KadB 

that is log-concave in z since ddB is positive . For maximal 
ratio combining of -branches i.i.d. Rayleigh fading, the PDF 
of at the combiner output is given by 

1 



/a=(e)= 



from which by 



(A^-1)! 
we obtain 

1 



h{z) = 



-,Nz 



(27) 



(28) 



(A^-1)! 

that is log-concave in z since N is positive and Nz — is 
concave. 

It is also important to remark that the log-concavity property 
for the EP proved above can have several applications, thus its 
relevance is beyond what illustrated in the following sections 
where an application example for bounds and local bounds is 
provided. 

III. Bounds and Local Bounds on Log-Concave 
Error Probability 

In this section it is shown how to take advantage of the 
log-concavity property. for the derivation of bounds and local 
bounds which are analytically simple and invertible for further 
analysis. An example of application will be briefly discussed, 
addressing local bounds of relevant performance metrics for 
adaptive A/-QAM systems. However, the application of the 
log-concavity is not limited to these cases (e.g., bounds for 
multidimensional modulation|3 as well as for A/-PSK can 
also be derived). 

The main idea to be exploited is that, due to the log-concave 
behavior proved in Sect. [Ill the EP plotted in logarithmic scale 
versus the signal-to-total disturbance-ratio, 7, in deciBel (dB) 
is a concave function (see, e.g., Fig|2]i. After having identified 
the ROI, where the system typically operates, we aim to easily 
obtain tighter analytically tractable and invertible upper and 
lower bounds valid in the ROI. The ROI is defined as the 
range [Pem, PeM\ of the EP which is of interest in the specific 
application. 

With the purpose to make a concrete example, in the 
following we consider the case of AWGN plus fading channel 
in which the performance is defined in terms of mean EP, the 

'^See , e.g., [33], [34]. The benefit provided by multidimensional constel- 
lations has been widely known in the design of coded modulation [35], [36]. 
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Fig. 2. A general behavior for the EP in log-scale versus 7 in dB: concave 
exact solution, upper and lower bounds (UB, LB) as well as local upper and 
lower bounds (LUB, LLB) in the ROI of interest are reported. 



EP hereafter, averaged over small-scale fading as a function of 
the mean 7, that is 7. Since the EP is monotonically decreasing 
in 7, the ROI corresponds to the SNR range [7Mi7m]' with 

Peir'J = ^em and P,{^1,) - P^M- 

Let us first consider bounds valid for all SNRs, that is 
for a ROI corresponding to SNR in the range [0,oo]. This 
ROI includes asymptotic behavior of EP. As an example, it 
is worthwhile to recall that in several cases, such as in single 
and multiple channels reception fading channel with Rayleigh 
or Nakagami-m PDF, the system achieves a diversity V if 
the asymptotic error probability is log-linear. This means that 
Pe(7) w ^(7) = K/^"^ where is a constant depending 
on the asymptotic behavior. In other words, a system with 
diversity T) is described by a curve of error probability with 
a slope approaching lO/P [dB/decade] for large 7. 

Thus, we focus our attention on systems with log-linear 
asymptotical mean EP. For these systems the log-concavity of 
the mean EP immediately implies that its asymptotic behavior 
provides an upper bound in the ROI [0, cx)]. Let us consider the 
usual EP, Peij), and the asymptote ■i;(7), both in logarithmic 
scale as a function of 7 (in dB). Note that on this scale the 
EP is concave whereas v{^) is linear It is clear that Pe(0) < 
v{Q), since the EP is less than or equal to 1/2. Furthermore, 
since Pe(7) and ^(7) are both decreasing, and the two curves 
approach at 7 — > 00, then Pe(7) < "^(7)- Therefore, an upper 
bound to Pe(7) can be easily defined as 



Pev^il) = min \ ^,w(7) 



(29) 



The UB on the inverse EP, that is on the value of 7 required 
to reach a target EP, P* < 1/2, is thus given by 

i/v 



K 



(30) 



To define local bounds let us now consider a generic ROI 
[Pem^Peu] with < Pem < PeM < 1/2. By shifting the 
asymptotic UB to touch the exact solution at extremes of the 
ROI, we define local bounds tighter than previously known 
bounds, easily invertible, and thus enabling further analysis. 



At the lower end of the ROI, that is for a target EP P* = P^m, 
we can define A7„j(dB) as the difference, in dB, between 
the required 7 for the asymptotic upper bound and the exact 
solution: 

A7„,(dB)=7uB(dB)-7*(dB). (31) 



Then, in linear scale: 



A7. 



7uB(-fem) 



7^ 



We now define the LUB in the ROI as: 

^'eLUB(7) - min{PeM, ^'eUB(7 ' ^7m)} : 



(32) 



(33) 



which is an invertible upper bound within the ROI. In fact, 
for a target EP P* in the ROI, the LUB on the required 7 
becomes: 



7lub 



%b{p:) 
A7™ 



7n 



~p* 



(34) 



Thus, to define the LUB one needs to know the exact required 
SNR at one point, namely, at the lower end of the ROI. 

Similarly, one can define the invertible LLB, which is a 
lower bound within the ROI, by shifting the UB of A7^j 
referred to PpM- This needs only the knowledge of the required 
SNR for the EP at the upper end of the ROI, 7m. The LLB is 
given by: 



-PeLLB(7) - max{Pem, PgUB (7 ' ^1 Al)} ■ 



The LLB on the required 7 results in 



— * _ 7UB(-Pe ) _ — * 

7llb — — 7m 



A7Af 



1/V 



(35) 



(36) 



A qualitative example of bounds and local bounds within 
the ROI is reported in Fig|2] At this point it is important to 
emphasize that, while the LUB is, within the ROI, a tighter 
bound than the asymptotic UB and still invertible, the LLB 
obtained by translation of the UB can be worse compared to 
known LB, but in the other hand, the LLB is easily invertible 
enabling further analysis^ 

Remark: the log-concavity property opens the way for 
the definition of other classes of bounds, such as based on 
tangent in the extremes of the ROI or based on saddle-point 
(steepest descent method). Local bounds here proposed have 
the advantage of being simple and analytically tractable for 
further analysis. 

We now discuss a possible application of local bounds on 
direct and inverse BEP to the evaluation of relevant perfor- 
mance metrics for adaptive Af-QAM systems. Let us consider 
as an example, A/-QAM with coherent detection and A^- 
branches MRC, whose exact BEP averaged over i.i.d. Rayleigh 
fading is given in [6] and its asymptotic upper bound is in 
the form P^ub (7) = min {l/2,TAr(M)/7^} where Tn{M) 
depends only on the constellation-size and the diversity order. 
From PfeuB (7) one can obtain the upper bound on the inverse 

'^Asymptotic expressions for the EP in the form of v{j) can be found 
using "systematic" approaches when exact EP expressions are not available 
or asymptotic expressions can not be easily deduced from well-know (but 
often complicated) expressions (see, e.g., [37], [38]). 



CONTI, PANCHENKO, SIDENKO, & TRALLI: LOG-CONCAVITY PROPERTY OF THE ERROR PROBABILITY WITH APPLICATION 



7 



BEP, which is the bound on the SNR required to achieve a 
target BEP equal to < 1/2, and from this invertible LUB 
and LLB through ( [34] i and ( [36] ), respectively. This enables 
the derivation of LLB and LUB on the error outage (EO), 
outage probability based on the EP [30], [39], which is an 
appropriate QoS measure for digital mobile radio when small- 
scale fading is superimposed on shadowing (typically modeled 
with log-normal distribution [1], [40])[3ln systems with slow 
adaptive modulation}^ for a given target BEP, P^, the spectral 
efficiency (SE) is a discrete RV with distribution that depends 
on the SNR thresholds and on how they are computed (i.e., 
on the BEP expression of the given system configuration). Let 
Mj and 7*dBj be the j— element from the set of possible 
constellation sizes and corresponding SNR threshold (in dB), 
respectively, to achieve a target BEP. Then, the mean SE 
results in 

j-i 

3=0 

+ Mj[l-F^J%B,j)] , (37) 

where Mk = log2 M/j and F-^^{-) is the cumulative distribu- 
tion function (CDF) of = 101og]^Q7. By substituting in 
( |37] | the required SNRs, 7^3 ^ with LUB, 7lub dB j' '^^ obtain 
a LLB on the mean SE allowing a conservative design of the 
communication system with different constellation-sizes and 
diversity orders. 

IV. Conclusions 

In this work we proved an important property of the error 
probability as a function of signal-to-noise-ratio in dB for 
AWGN channel as well as AWGN plus fading channels 
with single and multiple channels reception. In particular, we 
proved that the error probability is log-concave for a wide 
class of multidimensional modulation formats which include 
A/-QAM for two dimensions. This property can have several 
applications. As example, we exploited log-concavity to derive 
upper and lower bounds and to define local bounds that are 
tight in a given region of interest for the error probability. 
We also discussed an application of local bounds highlighting 
the possibility of easy computation for the inverse of EP 
formulas without loosing significant accuracy in the evaluation 
of figures of merit interesting in wireless communications. 
However, we believe that the relevance of log-concavity prop- 
erty goes beyond the example provided in the paper and may 
be exploited for other different purposes. 
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small-scale fading [4], [6]. 



Appendix 

Proof: [Proof of Lemma [T] Inequality ( fT4l i holds for h = 
1, hence, from now on we will assume that < ft. < 1. First, 
let us prove this inequality in the one-dimensional case d = 1. 
In this case 

H{h) = P„M_l+p^/, 

H' (h) = + Pi 
H" (h) = . 

Therefore, we need to prove that {\ — h) H' > \ — H or 

{l-h)(^-Po + P^>P,^^ + P,{l^h), 

which results in the exact equality. Let us now consider the 
case d > 1. We write the factor 1 — iJ in the second term in 
( fT4b as 

l-H{h)=Y^P,[l-{^-±^yh^^ 
= Y,Pk{l-Hk{h)) 

k 

and let us think of the left hand side of ( fT4] l as a homogeneous 
quadratic form in {Pk)o<k<d of the type 

E PkPmHk,m (h) > 0, 
0<k,m<d 

where Hk^m{h) is given by 

Hk,,n (h) = {l-h) H't,{h)Hl^{h) (38) 
- [l-Hu{h)][H'^{h)-{l-h)H':^{h)] . 

Lemma [T] then follows from Lemma |2l □ 
Proof: [Proof of Lemma |2l Let us start by recalling 
the following well known identities involving binomial co- 
efficients: 




As for notation, if L is a linear combination of {Pk)o<k<d we 
denote with the coefficient of Pk in L. By definition of 

Hk,m in ( l39] l to finish the proof, it is enough to show that for 

any < k,m < n 

(1 - h) {H'}, {H'}^ >{l- H}^ {H' -{1-h) H"}, . 

(42) 

Since 

{H},=g''-^h^=(^y \\ (43) 
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^fd-m\^2{h + h^ + ... + h'^+") - {hd + m) + {d-m- 2a)) 

{l-h)J2('^ ''^){'^h{l + 2h+... + {m + a- 1) + ;i™+"d) (49) 

r{l) [a I (c^("^ + a-l)-2(l + 2 + ... + (m + a-l))) 

= j (^c? (m + a - 1) - (m + a - 1) (m + a)^ 

= ^ P ~ ™ j ('(-d -m + dm + m^) + a (rf + 2to - 1) + a^) (50) 



we have and ( [39] l. the left hand side of ( |46] | can be rewritten as 



TO 

a / ^ — ' \ a 



- MTTTT) + ^^^^ I mTTTT) j ^^^^ 



{k"^ -k) +2hk{d-l) + h'^d{d-l) 



d{\~h)Y,(^ J^\{h + e + ... + (47) 



/l2 (/l+l)' 

By plugging (ED, dll and dS into m we obtain ^^"^^1^1"!^' *^ "^^^^ '^^"'^ ^^^^^ ® ^^^^^ ^° 

(1 - H) + TO) ^ 2. E f ' - i (1 + + . . . + z^™--) 



(fc - fc2) - 2hk (d - 1) + /^^ ((2d - 1) fc + 2d - d^) + . _/d-TO\ 

X (2 (/i + + . . . + - (/id + m)) 



ft.d + fc 

The derivative of the right-hand side with respect to k is equal 
to 

+ Mi2k-l) + h^d id - 1)) ^ ^ 

which is negative if A: > 1/2. Therefore, the right-hand side /d — m\ 

attains its maximum for fc or fc = 1. The difference for ( I (d — to — 2a) — (48) 

fc = and A: = 1 is a ^ ^ 

— — — —h > after some mathematical manipulations we obtain ( |49] l. 

Finally, comparing expansions for the left and right hand 
and, thus, the maximum is attained at fc ^ and we need to ^j^jg^ ^^at is and dHJ, respectively, (|46j becomes 
prove that 



(1 - h) {g}„, jhd + m) {2d -d')+ h^d^ yf'^-'^^dil + h 

- hd ■ vv « y ' 



'm-\-a — 2\ 



This is equivalent to > 2 ^ (^'^ (l + 2/z + . . . + (to + a - 1) h"'+"-^) . 

hd (2'^-"' -{h + if'"" /i") > (46) 

2d-m ^ (h^ Y^d-m j^rn ^ Combining all the coefficients for each power of h, the left 

Y~7h {h + lf^"'h"' (hd + m) . hand side can be written as r (h) = cq + cih + . . . + Cd-2h'^~'^ 

where 

Using the fact that 

d—m /J 

d — m\,„,. „ _ /-J o/'; I I a — m 



ri a—m 
(/l+l)'^-™/l"= ^ Q = (d-2(Z + 1)) 



Q>(2+i-m)V0 



a 
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Notice that the sign of c; is determined by d—2{l + l) so that 
c; > if Z < Zq and q < if Z > Iq. Since h < 1, this gives 

r{h) = co + cih + ... + Cd-2h'^^^ 
= {co + cih + . . . + cigh^") 

> {cQ + ci + ... + ci,)h'° 

+ + l + i + . . . + Cd-2) 

> h'" {C0 + C1 + ... + Cd-2) = h'r (1) . 

It remains to show that r (1) > 0. We observe that r (1) results 
in (|50l). Then, using ([39l), (gO]), and (STJ we obtain 

r (1) = 2'^-"'-^ [4 (-d - m + dm + m^) 

+ 2 (d + 2m - 1) (rf - to) + (d - TO - 1) (d - m)] 
= 3(i^ — 5(i + 3((i— l)m + dm + 

which is positive for d> 1. □ 
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